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One of the most important prerequisites for defining the spectrum of a 
nonstationary process is that the spectrum should transform simply and 
reasonably when the process is transformed linearly, and should lead to 
information about the important response statistics. Presented in this paper 
are some useful transform relationships for linear causal systems in terms 
of C. H. Page's time-varying spectra. Expressions suitable for direct 
analysis or numerical compulation of the time-history of the response 
process and its bounds, the response power spectrum, the total energy of the 
system, and the upper bounds on the response shock spectra are given. 

I. INTRODUCTION 

R. M. Loynes 1 recently established a list of desirable properties for 
the spectrum of nonstationary processes. The elementary properties 
include: (i) a nonstationary spectrum should be rigorously denned, 
(ii) it should describe in some sense the energy distribution over fre- 
quency and time, and (in) it reduces to the ordinary spectrum when 
the process is stationary (Loynes' properties Al, A2, and A5). Both 
Page's instantaneous power spectrum 2 and M. B. Priestley's evolu- 
tionary spectrum 3 satisfy these basic requirements. Another spectrum 
definition based on the notions of two-dimensional spectra arising in 
the consideration of harmonizable processes 4 - 5 also satisfy these basic 
requirements with some qualifications. However, from the practical 
point of view, especially when the filtering and convolution of a random 
process is involved the most important requirement is the existence 
of simple transform relationships for linear systems. A spectrum should 
transform simply and reasonably when the process is transformed 
linearly (property A3, Loynes). In other words, input-output rela- 
tionships are required so that a knowledge of the spectrum of a process 
determines the spectrum of the transformed process. In addition, these 
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relationships must be suitable for numerical computations. This is 
because most vibration data are available in the discrete form and 
the analytical solution of most nonstationary problems is difficult to 
obtain. Presented in this paper are some useful transform relationships 
for nonstationary processes using Page's time-varying spectra. 

II. BACKGROUND 

Consider a real function x(t) in (— oo, oo), which may be either a 
sample function of a random process, or a shock function which is 
zero outside the range t <. [a, b] for finite a and b and is Riemann 
integrable in this same range. The function x(t) is assumed to have 
finite energy. The running spectrum of x(t) is defined as 2 

X(t, to) = / x(t) exp ( — iior) dr, (1) 

J— 00 

where "i" is the imaginary unit. The instantaneous power spectrum 
is defined as 

Px (t, co) = f t | X(t, w ) | 2 

= 2x(t) Re [exp (iat)X(t, «)] (2) 

which is even in u>. If x(t) is a random process having a time- 
dependent autocorrelation function R x (t, t) = E[x(t)x(t — t)] where 
E denotes expectation, the instantaneous power spectrum of the process 
is understood to be the average of the spectra of all its sample func- 
tions. LetS x (t, <o) = E[p x {t, o)) ] and take the expectation of both sides 
of (2) , we obtain 

S t (t,to) = 2 [ R x (t, t) cos o>t dr. (3) 

The above relation shows that S x (t, w) is completely defined by 
R x (t, t). When x(t) is stationary, (3) reduces to the ordinary relation 
as _R.rU, t) = Rx(r) and S x (t, w) = S(a>), both independent of t. 

Some important properties of X(t, co) and p(t, w) as evident from their 
definitions are in order. 

(i) X(t, u) is Hermitian, i.e., X(t, -co) = X*{t, co). The symbol * 
denotes the complex conjugate. 

(ii) Let E{t) be the energy of x(t) up to time t. The function p x (t, w) 
can be regarded as the energy density in the (t, co) plane as seen by 
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E(t) - f \x(r)\ 2 dr = f ^ | X(t, «)| a &, 

= £ /^ /_" Px(r, c) rfco dr. (4) 

(m) Both X(7, w) and p x (<, w) depend on the past history of x(t), 
but not on the future. 

(iv) The function p z (t, u>) may take negative values (see Ref. 2, p. 106) 
but the integral /!_„ p z (r, 10) dr, the energy spectral density of frequency 
co, is always positive. 

(v) The time derivative of X(t, w) relates to the signal x(t) itself by 

x(t) = exp (itat) dX(t, u)/dt 

= exp (-uat) dX*(t, <a)/dt. (5) 

We will now proceed to derive from some input-output relations for 
a simple causal linear system in terms of X(t, <■>) and p x (t, w). 

III. INPUT-OUTPUT RELATIONS 

Consider a simple, second-order causal linear system whose equation 
of motion is given by 

y + 2\u y + uly = x(t) (6) 

with homogeneous initial conditions y(Q) = y(Q) = 0, where X < 1.0 
and u are positive constants representing respectively the damping and 
natural frequency of the system. This system has a transfer function 
between the output y and input x given by H(co) = (u 2 — to 2 + 2z'Xww ) -1 
and a corresponding impulse response h(t) = [exp (— \io t) sin pt]/p for 
t > 0, where p = (1 — X 2 )^ . We obtain the following results for 
different response parameters. 

3.1 Time-Varying Spectra oj Response Process 

Using the convolution integral and the casual property of the sys- 
tem, i.e., h(t) = for t < 0, and assuming the excitation begins at 
t = 0, the response y(t) in (6) is given in terms of X(t, u) by 

y(t) = [ x(r)h(t - r) dr 
Jo 

= y f H(cS)X(t, to) exp (uot) die. (7) 

Following the definition of (1), the running spectrum of y(t) is 
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Y(t, co) = X(t, co)#(co) - e Y (t, co), (8) 

where 

e Y (t, co) = / x(t,) exp (— iwTi) I h(r 2 ) exp (— iur 2 ) dr 2 dry (9) 

•'O Jt-Tl 

is the error involved in the approximation Y(t, <o) a: X(t, <a)H{o>). 
Notice from (8) and (9) that Y(t, o) is given by the double integral 
in the n — t 2 plane of a function /(n , t 2 ) = xMhfa) exp (— iuti 
— iwt 2 ) over an area of an equal-sided triangle formed by the line 
Tl + T2 = t in the first quadrant. The value of X(t, a>)H(u>) is given 
by the same integral over an extended area from to t along the 
ti axis and from to oo along the t 2 axis. The error e Y is given by the 
difference of these two integrals. 

Let F{x t ) be the Fourier transform of x t , the part of x(t) from 
t on, and F{y t ) be the Fourier transform of y t , the part of y(t) from 
t on. It can be shown that the error e Y {t, <a) in (9) is also given by 

e Y (t,co) = H(o)F[x t \ -F\ yi }. (10) 

It follows from (8) that the instantaneous power spectrum of y(t) 

is 

*(«. -) = (f B - %>*•** - •» + <** - 4^ ** - f ) 

= p,(«,«) |H(«) !■ + «,(«,«), (11) 



where 



,, ((l-) .l|. rl .-ir e ,f _i^ r ^ (12) 



is the error involved in the approximation p u (t, co) = p z (t, co) | i/(co) | 2 . 
Some interesting remarks can now be made about e Y , and subse- 
quently about e„ which is closely related to e Y through (12). The energy 
bound in co for e Y (t, co) can be considered as follows. Let the norm of an 
arbitrary function Q(co) be defined as || Q || u = [/"„ | Q(co) | 2 dco]*, then 
from Parseval's Theorem || F{x t } || u = || x t \\ T = [/" | x(t) \ 2 d t]*, and 
similarly || F\y t ) \\ m = \\ y t || T = [/" \y(r) | 2 drf. Suppose h(t) is inte- 
grable so that | //(co) | S A is bounded, then it follows from (10) that 



er(t, co) |L ^ A 



/^ oo ~u r /.oo "i| 

| X(r) I 2 dr J + j t | y(r) | 2 rfr J 



which — > as £ — > °o. If furthermore ,r(<) is integrable, from which 
| F\x t ) | ^ J1 | z(t) I dr which -> with i and consequently | F{?/,} | 
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also — » with t, then (10) indicates | e y (t, o>) | — > in t, uniformly in to. 
This result can also be concluded from (9) because in which the function 
/(ti , t 2 ) = x(t 1 )Ii(t 2 ) exp ( — koti — tor 2 ) is integrable in two dimensions. 

For a more practical purpose, we now consider the familiar first- and 
second-order simple linear systems. The running spectrum of the 
response y(t) and other pertinent quantities are given in Table I, in 
which a -1 = t = decay time of the impulse response h(t) of the first- 
order system and t = (Aw ) _1 decay time of h(t) of the second-order 
system. From this table, it can be seen that for large t/r ratio the 
error ey(t, w) for both systems [given by the second term in the ex- 
pression for Y(t, a))] will be small. In general at the same time t, the 
error for a system having short decay time (broad bandwidth) is 
smaller than that having long decay time (narrow bandwidth) . There- 
fore e Y (t, w) for systems with a flat spectrum (high damping) is 
smaller compared with that for high resonant systems (low damp- 
ing) . This observation is also evident from the damping term appear- 
ing in the expression of Y (t, <o) for the second-order system in Table I. 

In many cases the impulse responses and the input are both 
bounded functions defined on positive t axis such that \h(t)\ 52 
Ai exp(— ait) and \x(t)\ 2| A 2 exp(— oat) for some positive At , Az , 
and a x > a 2 . Then from (9) \e Y {t, «i) | ^ AiA 2 [exp {—a 2 t) — 
exp {—a l t)]/a 1 {ai — a 2 ). Therefore the error \e Y {t, w) | approaches 
zero as the time t increases ; and for a certain time t > t n the error can 
be regarded as negligibly small. 

From (8) another useful pair of equations relating to the time 
derivatives of the input and output running spectra can be obtained: 



dY(t,u>) dX(t, u ) m , de Y (t, u ) 
dY*(f,u) dX*(t,o>) m( . de* Y (t, JO) 



(13) 



which in turn will lead to the following relationship between the input 
and output time histories. 

3.2 Response Time-History and Some Bounds 

Using (13) and the analogous relation for y(t) as x(t) in (5), we 
obtain 

y(t) =x(t)H(u)-^e i " t 

= x(t)H*(<*)-^re—' (14) 
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which relates the input and output time histories in terms of the 
transfer function of the system and the error function in (9). The 
time-history of the response function y(t) can also be related to the 
time-history of the input function x(t) through the simple expression 

x(t)y(t) = - f Re [H(u>)Mt, «) d». (15) 

IT J n 

Notice that the above relation is advantageous when the direct multi- 
plication of input and output functions is involved. We prove (15) in 
the following. Consider the integral 

/ x(t)v(t) dr = — I //(«) / x(t)X(t, u) exp (tW) dr du>. 

It follows from the causal property of the system and Parseval's 
theorem that 

Jf x(r)y(r) dr = ^ £ tf(o>) | X(t, W ) | 2 da. 

Differentiating both sides with respect to t, (15) immediately follows. 
We now establish the upper bound for \y(t) | and \y(t) | 2 in terms of 
transfer characteristics of the system and input spectrum. From (7) , 
Schwarz's inequality yields 

| y(t) | 2 g E(t)N, (16) 

| y(t) | ^ E(tfN*, (17) 

where E(t) is the total energy of the input function x(t) up to the time t 
and is given by (4), and N = j*„ \ h(r) | 2 dr > 0. Equations (16) and 
(17) are also valid if x(t) is a random process. In this case the left-hand 
side quantities \ y(t) \ 2 and \ y(t) \ are replaced by E[ | y(t) | 2 ] and 
E\ 1 2/(0 I ] respectively. For our simple system, N = (4XwJ) -1 , and 

^ [| y(t) |2] - g^xs /_! /_" Si(t ' w) dw dr - (i8) 

3.3 Total Energy of the System 

The total energy E s {t) of the system described by (6) at time t is 
given by 

E,(t) = \ | Z{t, p) | 2 exp (2Xw„0. (19) 

where Z(t,p) is the running spectrum of «(i) = as(£) exp(Aw t) at the 
frequency p. The relation (19) can be proved as follows. Taking, 
according to (1), the running Fourier transform of both sides of (6) 
and integrating by parts, we have 
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[y + O'o> + 2\u )y] exp (-iut) + [<a\ + 2z\coco - u 2 ]Y(t, a>) = X(t, cj) 

in which the initial conditions have been used. The second term on the 
left-hand side of the above equation vanishes when setting w s $ ± = 
i\o> + p or a) = 9-2 = iAwo — p. Thus 

[0 + (»0i + 2Xco )2/] exp (-*^0 = X(«, 50 = Z(«, p), (20) 

[2/ + (tfi 2 + 2Xw.)y] exp (-iM - X(«, 2 ) = £*(*, p). (21) 
Multiplying (20) and (21) and dividing the result by 2, we have 

My 2 + 2X«.y0 -I- «y] = 1 1 Z(t, p) | 2 exp (2Xo>.«). 
But the left-hand side is the total energy of the system and therefore 
(19) is proved. 

Taking expectation of both sides of (19), we obtain the expected 
energy of the system subjected to the random excitation x(t) 



E[EM = \j j R*(r* . ^ ~ *•*) 



• exp [Xco„(t, + t 2 ) — ipfa — t 2 )] dr x dr 2 . (22) 

For the undamped system, i.e., at the limit when X — » 0, (19) reduces 
to \{y 2 + 2/ 2 ) = \ I 3T(<> «o)-| a which agrees with the principle of energy 
conservation. It is also interesting to note that a simple manipulation 
of (20) and (21) can lead to the familiar convolution relation between 
input and output 

1 f' 
2/(0 = " / X ( T ) ex P [-H0 _ T )] smp(t — t) dr 

P J -00 

as given by (7). 

3.4 Bounds on the Shock Spectrum 

The shock spectrum of x{t), 1 1 [a, b] is the maximum absolute response 
defined as S d (u , X) = sup, \y\ = the displacement spectrum or 
S,(o} , X) = SU P< \il\ = * ne velocity spectrum. We establish the 
following upper bounds to them: 

&(«„ , X) ^ sup - | Z(t, p) | exp (-Xw o U sup - | Z(t, p) | , (23) 
t LP J i P 

S v (<»o , X) ^ sup L ^ x y | Z(t, p) | exp (-Xw.Oj ^ sup | Z(«, p) |. 

(24) 
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The proof of these relations is straightforward. Write Z(t, o) = 
\Z(t, a>)\ exp [— i(j> z (t, to)] in which the phase angle <j>~(t, w) = tan -1 
[— S~(t, <a)/C g (t, <o)],S B and C s are the running sine and cosine trans- 
form of z{t) respectively. It follows from (20) or (21) that 

y + \u„y + ipy = | Z(t, p) | exp {-\u>J + i[pt - <t>,(t, p)]\- 
The above relation leads to 

y + Xoj,,// = | Z(t, p) | exp (-Xa> o cos [pt - <t> z (l, p)], 

VI) = I Zfo V) I exp (-Xw„0 sin [pt - <t> z (t, p)], 

and (23) and (24) follow directly. 

For the undamped system (A — > 0), p = o> , Z(t, p) = X(t, m ) ; 
therefore (23) and (24) reduce to 

a) a S,,(u , 0) ^ sup | X(t, wo) | and S,(w„ , 0) ^ sup | X(t, w ) |, 



i 



respectively. 



IV. ILLUSTRATIONS 



We now consider the following nonstationary processes for their 
time-varying spectra. 

4.1 Multiplicative Process 

x(i) = n(t)4>(t), 

where n(t) is a stationary process with an autocorrelation function 
R„(t) and spectral density S„U>), and <f>(t) is a causal deterministic 
function. The time dependent correlation function and power spectrum 
are 

R x (t, t) = R n (r)<t>(t - r), 

S t (t,<a) = 2<t>(t) [ R n (T)<p(t - t) coscorf/r. 
Jq 

When n(t) is a white noise with R„{t) = -R S(t), R > 0, then 
S x (t, w) = S(t) = 2R <f> 2 (t) becomes frequency independent. When 
tf)(t) is slowly varying such that <f>(t — r/2) ~ <j>(t + t/2), then 
R T {t, t) = <f> 2 {t)R„( T ) and S s (t, <■>) = <t> 2 (t)S n (<»)- This implies that 
x{t) is locally stationary and its power spectrum changes with respect 
to time not in the general shape but in its area covered with the 
w-axis only. 
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4.2 Periodic Nonstationary Process 6,7 

In this class of processes, the time dependent autocorrelation func- 
tion is given as R x (t, t) = 2^ fc (-r) exp (2-irikt/T) for a constant T. It 
can be readily shown that 

SM, co) = ^ 2>*(o>) exp {2*ikt/T), 

where ¥ fc (<o) is the Fourier transform of ^Mt) and has the property 
that * (co) ^ and ¥ fc U) = y k {2-rrk/T — o>). Moreover, it is easy to 
show that the time-varying power spectral density of the response 
process y{t) is also periodic with the same period T: 

S„(t, «) = -f S &(«) exp (2Tikt/T) 

in which q k M = * k {<»)H{<»)H{2Trk/T - co). 

The application of the rest of the input-output relations of Section 
III to these two types of nonstationary processes for either direct 
analysis or numerical computation is a straightforward exercise. 
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